Spin-polarized transport through a single-level quantum dot in the Kondo regime 
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Nonequilibrium electronic transport through a quantum dot coupled to ferromagnetic leads (elec- 
trodes) is studied theoretically by the nonequilibrium Green function technique. The system is 
described by the Anderson model with arbitrary correlation parameter U. Exchange interaction 
between the dot and ferromagnetic electrodes is taken into account via an effective molecular field. 
The following situations are analyzed numerically: (i) the dot is symmetrically coupled to two fer- 
romagnetic leads, (ii) one of the two ferromagnetic leads is half-metallic with almost total spin 
polarization of electron states at the Fermi level, and (iii) one of the two electrodes is nonmagnetic 
whereas the other one is ferromagnetic. Generally, the Kondo peak in the density of states (DOS) 
becomes spin-split when the total exchange field acting on the dot is nonzero. The spin-splitting of 
the Kondo peak in DOS leads to splitting and suppression of the corresponding zero bias anomaly 
in the differential conductance. 
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I. INTRODUCTION 



The Kondo phenomenon in electronic transport 
through artificial quantum dots (QDs) or single molecules 
attached to nonmagnetic leads was predicted theoret- 
ically more than a decade ago 0. Owing to recent 
progress in nanotechnology, the phenomenon has been 
also observed experimentally 0, Q • Several theoretical 
techniques have been developed to describe this effect 
SEIlllllHiniini- The description is usually 
simpler in the linear response regime, where equilibrium 
methods can be applied, but it becomes more complex 
when the system is driven out of equilibrium by an ex- 
ternal bias voltage 0, El El El El El One of the 
methods used to describe non-equilibrium Kondo effect is 
the non-equilibrium Green function technique d E, . 
To calculate density of states (DOS) and electric cur- 
rent one then needs the retarded/advanced as well as the 
lesser (correlation) Green functions. These can be de- 
rived within some approximation schemes. 

It has been shown recently that the Kondo effect can 
also occur when repla cing nonmagnetic leads by ferro- 
magnetic ones [la Gil l2il L~f- but ferromagnetism 
of the electrodes generally suppresses the effect - either 
partially or totally 0,Hj|. However, in some peculiar 
situations the effect remains almost unchanged. Sup- 
pression of the Kondo anomaly is a consequence of an 
effective exchange field due to coupling between the dot 
and ferromagnetic electrodes. The exchange field gives 
rise to spin-splitting of the equilibrium Kondo peak in 
DOS, and the two components of the splitted peak move 
away from the Fermi level, which leads to suppression 
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of the Kondo anomaly in electrical conductance - simi- 
larly as an external magnetic field suppresses the effect 
in nonmagnetic systems. Such a suppression was studied 
recently by the Green function technique in the limit of 
infinite correlation parameter U |19|. and was also con- 
firmed by numerical renormalization group calculations 
[23L |24|. However, only QDs symmetrically coupled to 
two magnetic leads were studied up to now. The Kondo 
anomaly survives then in the antiparallel magnetic con- 
figuration and is significantly suppressed in the parallel 
one. Recent experimental observations on Cqq molecules 
attached to ferromagnetic (Ni) electrodes support these 
general theoretical predictions |25j . 

Some features of the non-equilibrium Kondo phe- 
nomenon in QDs coupled to ferromagnetic leads have not 
been addressed yet. Therefore, in this paper we consider 
a more general situation. First of all, we consider the case 
when the two ferromagnetic electrodes are generally dif- 
ferent. In other words, the dot is (spin-) asymmetrically 
coupled to the ferromagnetic leads. This leads to qual- 
itatively new results. Second, we consider the case of 
arbitrary U instead of the limiting situation of infinite U 
studied in [2] . Third, we introduce an effective exchange 
field to describe the dot level renormalization. 

We analyze in detail three different situations. In the 
first case the dot is coupled to two ferromagnetic leads, 
and the coupling is fully symmetric in the parallel mag- 
netic configuration. We show that the equilibrium Kondo 
peak in DOS is then spin-split in the parallel configura- 
tion, whereas no splitting appears in the antiparallel one. 
The splitting, however, is significantly reduced for small 
values of the correlation parameter U. The correspond- 
ing zero-bias anomaly in conductance becomes split in 
the parallel configuration as well El- The second situ- 
ation studied in this paper is the one with asymmetric 
coupling to two ferromagnetic leads. As a particular case 
we consider the situation when one of the ferromagnetic 
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electrodes is half metallic, with almost total spin polar- 
ization of electron states at the Fermi level. Such struc- 
tures have been shown recently to have transport char- 
acteristics with typical diode-like behavior pa, [271 • Fi- 
nally, we also analyze the case when one of the electrodes 
is nonmagnetic whereas the second one is ferromagnetic, 
and show that one ferromagnetic electrode is sufficient to 
generate spin-splitting of the Kondo anomaly. 

The paper is organized as follows. The model and 
method are briefly described in Sections 2 and 3, respec- 
tively. Numerical results for the three different situations 
mentioned above are presented and discussed in Section 
4. Summary and general conclusions are given in Section 
5. 



II. MODEL 

We consider a single-level QD coupled to ferromagnetic 
metallic leads (electrodes) by tunnelling barriers. We 
restrict our considerations to collinear (parallel and an- 
tiparallel) magnetic configurations and assume that the 
axis z (spin quantization axis) points in the direction of 
the net spin of the left electrode (opposite to the cor- 
responding magnetic moment). Antiparallel alignment 
is obtained by reversing magnetic moment of the right 
electrode. The whole system is then described by Hamil- 
tonian of the general form 

H = H L + H R + H D + H T . (1) 

The terms Hp describe here the left (/3 = L) and right 
(/3 = R) electrodes in the non-interacting quasi-particle 
approximation, Hp = Y.ka e ^^ c %ia c ^!^^ wner e s k p a is 
the single-electron energy in the electrode /3 for the wave- 
number k and spin a (a =f , j), whereas ctg CT and c k p a are 
the corresponding creation and annihilation operators. 
The single-particle energy e k p a includes the electrostatic 
energy due to applied voltage, e k p a = e° kfj(T + eU{? = 
£ k/3<7 where s k p a is the corresponding energy in the 

unbiased system, U£ is the electrostatic potential of the 
/3-th electrode, e stands for the electron charge (e < 0), 
and up is the chemical potential of the /3-th electrode 
(the energy is measured from the Fermi level of unbiased 
system). Electron spin projection on the global quan- 
tization axis is denoted as | f° r s z — 1/2 and J, and 
for s z = —1/2. On the other hand, spin projection on 
the local quantization axis (local spin polarization in the 
ferromagnetic material) will be denoted as + for spin- 
majority and — for spin-minority electrons, respectively. 
When local and global quantization axes coincide, then 
t is equivalent to + and | is equivalent to — . (Note, that 
the local quantization axis in the ferromagnet is opposite 
to the local magnetization.) 

The term Hp in Eq.(l) describes the quantum dot and 
takes the form 

H D = ^T / e*d+d <J + Ud+d^d+d l , (2) 



where e CT denotes energy of the dot level (spin-dependent 
in a general case) , U denotes the electron correlation pa- 
rameter, whereas d^ and d a are the creation and anni- 
hilation operators for electrons on the dot. The level en- 
ergy e a includes the electrostatic energy due to applied 
voltage, e a = e 0a + eUf, where Uf is the electrostatic 
potential of the dot, and eo a is the level energy at zero 
bias. 

The electrostatic potential Uf of the dot will be de- 
termined fully self-consistently from the following capac- 
itance model [23, H3 : 

e (E n ° - E n 0°) = C L( U e - U e) + C R {U d e - Uf), 

F ( 3 ) 
where n a and no a are the dot occupation numbers {d+d a ) 
calculated for a given bias and for zero bias, respectively, 
whereas Cl and Cr denote the capacitances of the left 
and right tunnel junctions. Self-consistent determina- 
tion of the dot electrostatic potential makes the descrip- 
tion gauge invariant. This is particularly important for 
strongly asymmetric systems. 

The last term, Ht, in Eq.(l) describes tunnelling pro- 
cesses between the dot and electrodes and is of the form 

h t = Y, V wApJ° + h - c - » ( 4 ) 

fc/3<T 

where V k p a arc the components of the tunnelling ma- 
trix, and h.c. stands for the Hermitian conjugate term. 
Hamiltonian (4) includes only spin-conserving tunnelling 
processes. 

III. THEORETICAL FORMULATION 

Electric current flowing from the /3-th lead to the quan- 
tum dot in a nonequilibrium situation is determined by 
the retarded (advanced) and correlation (lesser) G< 
Green functions of the dot (calculated in the presence of 
coupling to the electrodes), and is given by the formula 

M 

^ = jj ^(E){G<(E) + fp(E)lG:(E)~G a a (E)}}, 

(5) 

where fp{E) is the Fermi distribution function for the 
/3-th electrode. The retarded (advanced) Green func- 
tions can be calculated from the corresponding equation 
of motion. The key difficulty is with calculating the lesser 
Green function G< ( E). 

In a recent paper |31| we applied the equation of mo- 
tion method to derive both G r a (E) and G<(E) Green 
functions within the same approximation scheme |32| . 
However, the approximations for the lesser Green func- 
tion G< (E) do not conserve charge current in asymmet- 
rical systems. Therefore, in the following we assume 
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constant (independent of energy) coupling parameters, 
T?{E) = 27rji Vk^V^J(E - e k(ia ) = T». As pointed 
out in Refs |3JJ, l34| . it is then sufficient to determine 
J (dE j '27i") G< (E) , while knowledge of the exact form of 
Gp(E) is not necessary. Current conservation condition 
allows then to express the above integral by an integral 
including retarded and advanced Green functions only, 
which in turn allows to rewrite the current formula in 
the commonly used form, 

^ = |/ ^^^{Gl(E)-G a a (E)}[f L (E)~f R (E)}. 

' 1 " (6) 
Similarly, the occupation numbers, n a = (d+d CT ), are 
then given by the formula 



dEY L J L {E)+Y*f R {E) 



2tt 



[Gl{E)-Gl{E)\. (7) 



In the following the parameters rg will be used to pa- 
rameterize strength of the coupling between the dot and 
leads. It is convenient to introduce the spin polariza- 
tion factors pp defined as pp = (T p + - Vt) / {T^ + Yt) , 
where and are the coupling parameters for spin- 
majority and spin-minority electrons in the lead /3, re- 
spectively. Accordingly, one may write = (1 + pp)r@ 

and r£ = (i - pp)rP, with r' 3 = (r£ + ri)/2. 

The retarded (advanced) Green function Ga^ of the 
dot can be calculated only approximately, for instance by 
the equation of motion method. In the approximations 
introduced by Meir et al @ one finds 



G r a (E) = 



E-e a -U{l-n_ a ) 



[E-e a - K(E)](E -e a -U)- Un_ a Z r a (E) ' 
where SJJ. is the corresponding self energy, 



*AE) - ^(E)+U LQa{E)[LQrT{E) _^ {E)] ■ 

(9) 

with Loo- = E — e a — U(l — n- a ), and T,q 1(7 (E) and 
Eq 3cf (E) defined as 

T,l 1<T {E)=n^l a {E)+Y,\ a {E), (10) 



^{E)=Y.UE) + Y,l a {E). 



(11) 



The self energies YT 0a {E), TT l(T (E), and T, r 3(J (E) are de- 
fined as 



i 



0=L,R 



0=L,R k 



E - e k p a + i0+ 



E 

0=LR 



de rg rg 

2tt E - e 2 



(12) 



and 



= E / d i A ^ 



P=lr: 



e — E — c — (j + e a — ih/i 



£ + E — e_ ff XJ ih/r- 



(13) 



(for i=l,3), where A\ = /(e), ^3 = 1, and t_ ct is the 
relaxation time of the intermediate states [j| . This relax- 
ation time is spin dependent and in the low-temperature 
limit is given by the formula [f| 



T = ^EE«:e (w 

Ta Z7Tn 0,0' <T< 



V0 



H>0' — Hp + £0 — e<r' 
(M/3 - - e<7') ! 



(14) 



where 0(x) = for a; < and Q(x) = 1 otherwise. 

In the limit of infinite {J the Green function (8) reduces 
to the well known form, 



Gl(E) = 



1 - n- a 



E-e a -^ a (E)-^ la {EY 
where Y> r 0rT {E) is given by Eq.(12) and Tl\ a (E) by 

£U£)= E&(25) 



(15) 



E 

0=L,R 



ds 



f0{eW- 



2-k-e + E- 



ih/i 



(16) 



and r CT defined by Eq.(14). 

The above derived Green functions are sufficient to 
describe qualitatively basic features of the Kondo phe- 
nomenon in QDs attached to nonmagnetic leads 
However, they are not sufficient to describe properly the 
Kondo phenomenon when the quantum dot is attached 
to ferromagnetic leads. The key feature of the system, 
which is not sufficiently taken into account is the split- 
ting of the dot level due to spin dependent tunneling 
processes |19| . The most natural way would rely on an 
extension of the Green function calculations by going be- 
yond the approximations used to derive Eq.(8). This, 
however, leads to cumbersome expressions. To avoid this, 
one may treat the problem approximately by introducing 
'by hand' the level splitting to the formalism described 
above. One way to achieve such an objective was pro- 
posed in Ref.[ljj for the limit of infinite U, where the 
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bare dot level e a in E^ cr (£') and T, r Q(J (E) was replaced by 
the renormalized energy e a calculated in a self-consistent 
way. Such a renormalization works well in the limit of 
infinite U . It is however not clear how to extend it to 
the general case of arbitrary U. Therefore, we decided 
to include the level splitting via an effective exchange 
field. The exchange-induced spin-splitting of the level is 
particularly important for the self-energies Yi\ a (E) and 
Yi^ a {E), so we replace the bare energy levels in the self- 
energies E \ a (E) and Eg CT (E) by the corresponding renor- 
malized levels ? CT = t a ± g^sB cx /2, with the exchange 
field calculated from the formula 

B cx = — £ Befp-Me) 



e-f — E7 — i^/r-f 




-0.10 -0.05 0.00 0.05 0.10 

eV/ D 



! 



£ - q - *Vn £ - e i - u - ih/n 



(17) 



When the bare dot level and spin relaxation time are 
independent of the spin orientation, e-j- = ej_ = e and 
T| = T| = r, the formula (14) acquires the form |3flj 
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FIG. 1: Bias dependence of the differential conductance in 
the parallel (a) and antiparallel (b) configurations, and the 
corresponding TMR (c) for indicated four values of the lead 
polarization Pl = PR = p (p = corresponds to nonmagnetic 
leads so the corresponding TMR vanishes and is not shown in 
part (c)). The parameters assumed for numerical calculations 
are: kT/D = 0.001, Y L /D = Y R / D = 0.1, e 0T / D = e 0l /D = 



e /D = -0.35. 
0.33. 



U/D = 500, and (e 2 /C L )/D = (e 2 /C R )/D = 



T 



e — e — ih/r £ — £ — [/ — ih/r 



(18) 

To some extent such an approach is similar to that 
used in Ref . 0] , and both approaches give similar results 
in the limit of large U. This follows from the fact that 
the expression (17) for exchange field is basically the ex- 
pression for the self energy Y7 la (see Eq.(13) for i = 1 and 
E = e a ). However, the approach based on the exchange 
field allows to handle easily also the general case of finite 
U. 



is included only via the spin asymmetry of the coupling 
parameters and T^. Apart from this, in all numerical 
calculations we assumed kT/D = 0.001 and eo^/D = 
eoi/D - e /D = -0.35. 

For positive (negative) bias we assume the electrostatic 
potential of the left (right) electrode equal to zero. In 
other words, electrochemical potential of the drain elec- 
trode is assumed to be zero while of the source electrode 
is shifted up by \eV\. In the following the bias is described 
by the corresponding electrostatic energy eV. Note that 
positive eV corresponds to negative bias due to negative 
electron charge (e < 0). 



IV. NUMERICAL RESULTS 

Now, we apply the above described formalism to the 
Kondo problem in a QD coupled to ferromagnetic leads. 
In the numerical calculations described below the energy 
is measured in the units of D, where D = D/50 and D 
is the electron band width. For simplicity, the electron 
band in the leads is assumed to be independent of the 
spin orientation and extends from ~25D below the Fermi 
level (bottom band edge) up to 25D above the Fermi 
level (top band edge). The energy integrals will be cut 
off at E = ±25 D, i.e., will be limited to the electron 
band. Thus, the influence of ferromagnetic electrodes 



A. QD coupled to two similar ferromagnetic leads 

Consider first the situation when both electrodes are 
made of the same ferromagnetic metal, and the coupling 
of the dot to both leads is symmetrical (in the parallel 
configuration). For numerical calculations we assumed 
= — 0.12 for spin-majority electrons and 

T^/D = T^/D = 0.08 for spin-minority ones, which 
corresponds to the spin polarization factor pl = pr = 
p = 0.2, and T L /D = Y R JD = 0.1. 

It was shown in Ref.|l9j that ferromagnetism of the 
electrodes leads to spin splitting of the Kondo peak in 
the density of states (DOS) in the parallel configuration, 
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whereas no splitting occurs for the antiparallel orienta- 
tion. Such a behavior of DOS has a significant influence 
on the transport properties. First of all, the Kondo peak 
in DOS leads to zero bias anomaly in the differential con- 
ductance Gdiff = dl/dV. This anomaly is particularly 
interesting in the parallel configuration, where the spin 
splitting of the Kondo peak in DOS leads to splitting of 
the differential conductance, as shown in Fig. 1(a) for four 
different values of the electrode spin polarization factor 
p and for large U. For p = there is only a single peak 
centered at zero bias. When the polarization factor be- 
comes nonzero, the peak becomes split into two compo- 
nents located symmetrically on both sides of the original 
peak, with the corresponding intensities significantly sup- 
pressed. The splitting of the Kondo anomaly increases 
with increasing p. Moreover, height of the two compo- 
nents of the Kondo anomaly decreases with increasing p. 
On the other hand, in the antiparallel configuration there 
is no splitting of the Kondo peak in the density of states 
and consequently also no splitting of the Kondo anomaly 
in the differential conductance (see Fig. 1(b)). For all po- 
larization values, the anomaly is similar to that in the 
case of QDs coupled to nonmagnetic leads. However, 
intensity of the Kondo anomaly decreases with increas- 
ing polarization. Difference between conductance in the 
antiparallel and parallel configurations gives rise to the 
TMR effect which may be described quantitatively by 
the ratio (I p - I AP )// AP , where I p and 7 AP denote the 
current flowing through the system in the parallel and 
antiparallel configurations at the same bias, respectively. 
The associated TMR effect is displayed in Fig. 1(c). One 
finds negative values of the TMR ratio, which is a conse- 
quence of the spin-splitting of the the Kondo peak in the 
parallel configuration and absence of such a splitting for 
antiparallel alignment. It is worth to note, that in the 
absence of the Kondo anomaly the TMR effect would be 
positive. 



The Kondo anomaly in transport characteristics shown 
in Fig.l was calculated for the limit of large U. In Fig. 2 
we show similar characteristics as in Fig.l, but for differ- 
ent values of the correlation parameter U and a constant 
value of the polarization factor p. The splitting of the 
Kondo anomaly in the parallel configuration and inten- 
sity of the peaks (Fig. 2(a)) decrease with decreasing U. 
In the antiparallel configuration there is no splitting of 
the Kondo anomaly, but intensity of the Kondo peak de- 
creases with decreasing U. The associated TMR effect is 
shown in Fig. 2(c). The effect is negative in a certain bias 
range around the zero bias limit, but absolute magnitude 
of the effect becomes smaller for smaller values of U. For 
large bias there is a transition from negative to positive 
TMR with decreasing U. 
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eV/D 

FIG. 2: Bias dependence of differential conductance in the 
parallel (a) and antiparallel (b) configurations, and the corre- 
sponding TMR (c) calculated for indicated values of the cor- 
relation parameter U and for p — 0.2. The other parameters 
are as in Fig.l. 



B. QD coupled to one ferromagnetic and one 
half-metallic leads 

Assume now that one of the electrodes (say the left 
one) is made of a 3-d ferromagnetic metal, the second 
(right) one is half- metallic, and the total coupling to the 
latter electrode is much smaller than to the former one. 
This is reflected in the spin asymmetry of the bare cou- 
pling constants, for which we assume T+/D = 0.28 and 
Tt/D = 0.12 for the left electrode, and Tf/D = 0.04 
and T^/D = 0.0002 for the right one. These parame- 
ters correspond to p L = 0.4, p R = 0.99, T L /D = 0.2, 
and T R /D « 0.02. Thus, the spin asymmetry of the cou- 
pling to the right electrode is much larger than to the left 
one. In Fig. 3 we show DOS in the parallel (left column) 
and antiparallel (right column) magnetic configurations, 
calculated for vanishing as well as for positive and nega- 
tive bias voltages. Consider now the main features of the 
spectra in more details, and let us begin with the parallel 
configuration (left column in Fig. 3). 

At V = the Kondo peak in DOS is spin-split, and the 
intensity of spin-down peak is relatively large, whereas 
that of the spin-up peak is much smaller. The asymmetry 
in peak intensities is a consequence of the spin asymme- 
try in the coupling of the dot to metallic electrodes - this 
coupling is larger for spin-up electron, which determines 
hight of the Kondo peak for spin-down electrons. When a 
bias voltage is applied, each of the two Kondo peaks gen- 
erally becomes additionally split into two components. 
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FIG. 3: DOS for spin-up (thick lines) and spin-down (thin 
lines) electron states on the dot in the parallel (left column) 
and antiparallel (right column) magnetic configurations, cal- 
culated for three indicated voltages and for F+/D = 0.28, 
Vl/D = 0.12, Tf/D = 0.04, Y^/D = 0.0002, U/D = 500, 
(e 2 /Cl)/ D = (e 2 /Cr)/D = 10. The other parameters are as 
in Fig.l. 
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FIG. 4: Current-voltage characteristics in the parallel (a) and 
antiparallel (b) configurations, calculated for the parameters 
as in Fig. 3. Current in the spin-down (spin-up) channel in 
the parallel (antiparallel) configuration almost vanishes so the 
total current flows in the spin-up (spin-down) channel (the 
curves presenting the total and spin-up (spin-down) currents 
are not resolved. The insets show the corresponding occupa- 
tion numbers. The parameters as in Fig. 3 



One of them (the one associated with the coupling to the 
source electrode) moves up in energy, whereas position of 
the second one (the one associated with the drain elec- 
trode) remains unchanged. This is because we assumed 
that the electrochemical potential of the source electrode 
shifts up by |eV|, while of the drain electrode is inde- 
pendent of the voltage. For eV > (negative bias), the 
splitting of the large- intensity (spin-down) peak is clearly 
visible, although one component of the splitted peak is 
relatively small. This is just the component associated 
with the coupling of the dot to the right electrode in the 
spin-up channel. Since this coupling is relatively small, 
the corresponding intensity is also small. The second 
component, in turn, is much larger because it is associ- 
ated with the coupling to the left electrode in the spin-up 
channel, which is the dominant coupling in the system 
considered. Splitting of the low-intensity (spin-up) peak 
is not resolved. Intensity of the component associated 
with the coupling to the right electrode in the spin-down 
channel practically vanishes because this coupling is neg- 
ligible in the case considered. For eV < (positive bias), 
the situation is changed. Now the electrochemical po- 
tential of the left electrode is independent of the bias. 
Consequently, intensity of the components whose posi- 
tion is independent of energy is significantly larger than 
intensity of the other components (the ones associated 



with the right electrode). As before, the component as- 
sociated with the coupling to the right electrode in the 
spin-down channel is not resolved. 

Consider now the antiparallel configuration (right col- 
umn in Fig. 3), when magnetic moment of the right elec- 
trode is reversed. There is a nonzero spin splitting of 
the Kondo peak at equilibrium, contrary to the case of 
symmetric coupling to the magnetic electrodes, where 
the spin splitting in the antiparallel configuration van- 
ishes |l9j . Apart from this, the situation is qualitatively 
similar to the one for parallel configuration. The main 
difference is that now the bias-induced splitting of the 
large-intensity peak is not resolved, whereas splitting of 
the low-intensity peak is resolved. 

As in the case of symmetrical coupling described above, 
the Kondo peaks in DOS give rise to anomalous behav- 
ior of the corresponding transport characteristics. Due to 
the splitting of the equilibrium Kondo peak, the anomaly 
in DOS does not contribute to transport in the small 
bias regime. The Kondo peaks enter the 'tunnelling 
window' at a certain bias, which leads to an enhanced 
conductance. Such an enhancement is clearly visible 
in the current-voltage characteristics shown in Fig. 4 for 
both parallel (a) and antiparallel (b) configurations (solid 
lines), where for negative values of eV the enhancement 
is quite significant, but it is less pronounced for eV > 0. 
This asymmetry is due to the difference in intensities of 
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FIG. 5: Bias dependence of the differential conductance in 
the parallel (a) and antiparallel (b) configurations and the 
corresponding TMR (c), calculated for the same parameters 
as in Fig. 3. 
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FIG. 6: DOS for spin-up (solid lines) and spin-down (dot- 
ted lines) electron states of the dot, calculated for three 
different voltages, and for V%/D = 0.12, T*L/D = 0.08, 
Tf/D = F R /D = 0.1, U/D = 500, and (e 2 /C L )/D = 
(e 2 /Cr)/D = 0.33. The other parameters are as in Fig.l. 



the corresponding Kondo peaks that enter the 'tunnelling 
window'. 

The differential conductance in the Kondo regime is 
shown in Fig. 5 for parallel (a) and antiparallel (b) mag- 
netic configurations. In the parallel configuration the 
Kondo anomaly occurs in the spin-up channel and for 
eV > only. This may be easily understood by con- 
sidering the relevant DOS (see Fig. 3, left column). For 
eV > only the Kondo peak in spin-up DOS can enter 
the tunnelling window created by the bias. For eV < 0, 
on the other hand, the Kondo peak in spin-down DOS 
can enter the tunneling window. However, the spin down 
channel is almost non-conducting, so the corresponding 
peak in the differential conductance is suppressed. In 
the antiparallel configuration the Kondo peak in differ- 
ential conductance occurs for eV < only. This can 
be accounted for by taking into account behavior of the 
Kondo peaks in DOS shown in Fig. 3 (right column), and 
the fact that now the spin-up channel is non-conducting. 
For eV > only the Kondo peak in the spin-up DOS can 
enter the tunneling window, whereas for eV < this is 
the Kondo peak in spin-down DOS (of large intensity). 

The corresponding TMR is shown in Fig. 5(c). It is in- 
teresting to note that TMR is highly asymmetrical with 
respect to the bias reversal. It becomes positive for eV 
exceeding a certain positive value, and negative below 
this voltage. This is a consequence of the fact that for 
positive eV the Kondo peak in differential conductance is 
clearly visible in the parallel configuration (see Fig. 5 (a)), 
whereas for eV < the Kondo peak occurs in the antipar- 
allel configuration (see Fig. 5(b)). Such a behavior of the 



conductance and also TMR may be interesting from the 
point of view of applications in mesoscopic diodes. 

C. QD coupled to one ferromagnetic and one 
nonmagnetic leads 

A specific example of asymmetric systems is the case 
where one electrode is ferromagnetic (typical ferromag- 
netic 3d metal) whereas the second one is nonmagnetic. 
For numerical calculations we assumed T^/D = 0.12, 
T^/D = 0.08 for the left (magnetic) electrode, and 
= T^/D = 0.1 for the right (nonmagnetic) one, 
which corresponds to pl = 0.2, pa = 0, and T L / D = 
T R /D = 0.1. As in the other asymmetrical situations 
studied in this paper, the equilibrium Kondo peak in 
DOS becomes spin-split. When a bias voltage is applied, 
each component becomes additionally split, as shown in 
Fig. 6. Variation of the spectra with bias voltage can be 
accounted for in a similar way as in the case of the dot 
coupled asymmetrically to two ferromagnetic electrodes. 
The only difference is that now all components of the 
peaks are clearly resolved. This is because all coupling 
constants are now of comparable magnitude. 

The corresponding differential conductance is shown 
in Fig. 7. Due to the spin splitting of the Kondo peak in 
DOS, the Kondo anomaly in the conductance becomes 
split as well, as clearly seen in Fig. 7. However, the 
splitting is asymmetric with respect to the bias reversal. 
Thus, there is no need to have two ferromagnetic elec- 
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FIG. 7: Bias dependence of the differential conductance, cal- 
culated for the parameters as in Fig. 6. 



trodes to observe splitting of the Kondo anomaly, but it 
is sufficient when only one lead is ferromagnetic. 



V. SUMMARY AND CONCLUSIONS 

In this paper we considered the Kondo problem in 
quantum dots coupled symmetrically and asymmetrically 
to ferromagnetic leads. As an specific example of asym- 
metrical systems, we considered the case when one elec- 
trode is ferromagnetic, whereas the second one is non- 
magnetic. 

We showed that ferromagnetism of the leads gives rise 
to a splitting of the equilibrium Kondo peak in DOS for 
all asymmetrical situations. This generally takes place 
for both magnetic configurations when the two electrodes 
arc different. The splitting in both configurations also oc- 
curs when both magnetic electrodes are of the same mate- 
rial, but the corresponding coupling strengths to the dot 
are different. Indeed, such a splitting in parallel and also 
antiparallel configurations was recently observed experi- 
mentally pl| . When similar electrodes are symmetrically 
coupled to the dot, the splitting occurs only in the par- 
allel configuration. An interesting conclusion from the 
experimental point of view is that the splitting also oc- 
curs in the case when one electrode is nonmagnetic. 

The spin-splitting of DOS can lead to characteristic 
splitting of the zero bias anomaly in electrical conduc- 
tance. This in turn can lead to negative (inverse) tun- 
nel magnetoresistance effect. In highly asymmetrical sys- 
tems TMR can change sign when bias voltage is reversed. 
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